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This article presents a new factorization algorithm based on the implementation of exponential sums 
using optical interference and exploiting the spectrum of the light source. Such a goal is achievable 
with the use of two different kinds of optical interferometers with variable optical paths: a liquid 
crystal grating and a generalized symmetric Michelson interferometer. This algorithm allows, for 
the first time, to find, in a single run, all the factors of an arbitrary large number A''. 



The factorization of large numbers A'' is a lot more 
difficult than the reverse operation of multiplying large 
prime numbers. This difficulty is at the basis of encryp- 
tion systems [T] . A more recent approach to factorization, 
proposed by W. Schleich, exploits the periodic properties 
of truncated exponential sums [2J |3J 3] of order j: 



op, 



A 



N 



1 



M + 1 



M 

E 



exp 



-27ri mP 



(1) 



where M-|- 1 is the number of phase terms in the sum (Af 
is called the truncation parameter), N is the number to 
be factored, and j and I are positive integers, with j > 1 
and I < £ < y/N, respectively. For j = 2, the truncated 
exponential sum reduces to a truncated Gauss sum [5]. 
If ^ is a factor of N, all the terms interfere constructively 
and the truncated exponential sum assumes its maximum 
value, i.e. 1. On the other hand, if £ is not a factor of 
N, the truncated exponential sum assumes a value less 
than one, because of the destructive interference caused 
by the rapid oscillation of the phases terms of order j in 
Eq. ([ij. Obviously, the more terms involved in the sum 
(i.e. the larger the truncation parameter M), the better 
we can distinguish between factors and non factors. It 
turns out that one needs at least M ~ ^-{/N terms in 
order to discriminate factors from non factors [4 . 

Truncated Gauss sums have been reproduced using dif- 
ferent techniques [5j |6l |7l |8] . Unfortunately, all these past 
experimental realizations present two common problems. 
The first one is that there is only one knob (one phys- 
ical parameter) to vary the global ratio N/£: the ratio 
N/£ is known before the experiment is run. The second 
problem is that it is necessary to run the experiment for 
each possible trial factor £, to find out which ones are the 
factors. 

We wish to present a new approach, based on optical 
interference, which allows, for the first time, to obtain all 
the factors of a large number A'^ in a single run, for any 
value of N. Such an approach is generalizable to the re- 
production of exponential sums of any order j, with sub- 
sequent reduction of the number of resources compared 
to the past realizations. 

Moreover, in this procedure, we will not take into ac- 
count the term with m = in the truncated exponential 
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FIG. 1: Theoretical model of a generic M-path optical inter- 
ferometer, where opm = m-' Nu is the length of the mth optical 
path with m — 1,2, M. The optical paths are represented 
by arrows, which length increase quadratically, in the case of 
Gauss sums {j = 2). The incoming electromagnetic field is 
given by a polychromatic plane wave and the spectrum of the 
outgoing field, given by the interference of all the M optical 
paths, is measured by a spectrometer. 



sum ([T]), because it corresponds to a null phase, no mat- 
ter if the trial factor £ is a factor or not. 

Our procedure is based on two basic simple ideas. The 
first idea consists of the use of an Af-path optical inter- 
ferometer. Such an interferometer needs to be able to 
reproduce, for a definite wavelength A, M spatial modes 
with phase terms of order j and add them coherently in 
order to reproduce the truncated exponential sums. The 
second basic idea consists of exploiting the spectrum of 
the incoming light in order to reproduce all the possible 
trial factors in a single run for any value of N, measuring 
the intensity of the outgoing light as a function of the 
wavelength. 

We will now explain, in detail, the physical working 
principle behind the two main ideas stated above. 

In general, an Af-path optical interferometer (see Fig. 
[T|), interacting with an incoming polychromatic plane 
wave, allows the coherent superposition of M electro- 
magnetic phase terms (modes). The optical phase of the 
generic term, associated with the wavelength A and with 
the m*'' optical path, with m = 1, A/, is given by: 
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is the m optical path, rim is the index of refraction of 
the material in the path, and dm is the length of the path. 

The truncated exponential sum ([T]), not including the 
term m = 0, is reproduced by the coherent superposition 
of all the M phase terms in Eq. if the following 

conditions are satisfied: 

0„(A) = 2W^, (4) 

for each m = 1,2,..., M. 

In order to achieve this goal, we use a polychromatic 
source so that we can exploit all the spectrum of the in- 
coming light in order to reproduce all the possible trial 
factors at the same time, through the following corre- 
spondence: 

A = iu, (5) 

where u is an appropriate unit of measurement of length. 
Such a unit is chosen so that all the possible trial factors 
can be experimentally reproduced exploiting all the range 
of Fourier modes A emitted by a polychromatic source 
S, with associated intensity |£'5(A)p large enough to be 
measured. The other important step consists in varying 
the M optical paths in Eq. ^ in order to satisfy the 
following conditions: 

opm = m^Nu, (6) 

for each m— 1,2,..., Ad. 

If both Eq. ([5| and Eq. ^ are satisfied, the magni- 
tude of the outgoing electromagnetic field is given by the 
superposition of all the M exponential phases in Eq. Q , 
for each of the Fourier modes A = £m in the spectrum of 
the source. 

At this point, we can use a spectrometer in order to 
measure the intensity of the outgoing electric field 

\E^{M^N,j)\^^\Es{\)A^^'''\\)\\ (7) 

as a function of the wavelength A = tu; where Es{\) is 
the amplitude associate to the Fourier mode A emitted by 
the source and ^^^^'■'''(A) is given by Eq. 1 with I = ^. 
Knowing the spectrum of emission |i?s(A)|^f the source, 
it is easy to extrapolate the modulo squared |^j^^'"'''(A)p 
of the truncated exponential sum versus the wavelength 
A = £u. The factors correspond to the wavelengths which 
allow constructive interference and so maxima in the ex- 
trapolated intensity spectrum. 

It is important to point out that our factorization ap- 
proach allows for the first time to achieve two important 
goals. The first one is the implementation of two indepen- 
dent parameters corresponding to N and £, respectively. 
Moreover this approach allows us to determine the fac- 
tors of TV in a single run of the experiment. For the first 
time, in an exponential sum approach to factorization, 
the ratio y is not known a priori and it is not necessary 
to run the experiment for all the possible trial factors. 



At this point, it is important to point out some exper- 
imental aspects associated to the conditions in Eq. ([s]) 
and Eq. ^ . Both the optical paths and the wavelengths 
in the spectrum of the outgoing field need to be measured 
with the resolution of lu. So, experimentally speaking, 
the larger is the unity of measurement u, the smaller is 
the necessary resolution in the actual experimental real- 
ization. But, at the same time, increasing the value of 
u, the bandwidth of the optical spectrum necessary to 
cover all the trial factors increases by the same factor. 
Another important point to take into account is that the 
less the resolution of the spectrometer, the larger is the 
number of measurements we need to perform in the dif- 
ferent ranges of the spectrum of the light source, in order 
to cover all the trial factors. 

Now we present a simple factorization algorithm, based 
on the factorization approach described so far. Such an 
algorithm can be implemented in four different steps. 
The first one consists of choosing a particular wavelength 
A, which corresponds to a given number iV = ^ to factor- 
ize. The second step consists of varying the optical paths 
in Eq. ([s]) as long as the electromagnetic phase terms in 
Eq. ([2]) satisfy the following condition: 

</)„,(A)=2W, (8) 

with m = 1,2, ...,M. 

In the third step, we consider a polychromatic plane 
wave interacting with the generic M-path interferometer 
in Fig. [1} keeping fixed the optical paths determined 
in the previous step. In this case, each wavelength A 
corresponds to a particular set of phases Q , which can 
be written in the following way: 

^mW = 2Trm^ J, (9) 

with m = 1,2,...,M. The phases ^ are the exponen- 
tial sum phases, associated with the number iV = A to 
factorize and the generic trial factor £, given by Eq. (|5|. 

In the fourth step of the algorithm, we implement the 
factorization of a generic number N, choosing a parame- 
ter a such that: 

N ^aN = a-. (10) 
u 

In order to achieve this goal, it is possible to rescale the 
condition A = £u, imposing the new condition X = £^. 
Unfortunately, this approach is limited by the achievable 
range of wavelengths in the spectrum of emission of the 
source. Another possible approach, independent from the 
spectrum of emission of the source, consists of multiply- 
ing the lengths of all the M optical paths, found in the 
second step, by a factor a. In this way, we can factorize 
any number N using a simple rescaling procedure. 

Let us now describe how to experimentally realize our 
factorization procedure. In order to satisfy the condi- 
tion ([g]), we need to be able to manipulate either the 
indexes of refraction associated with the M optical 
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FIG. 2: Voltage-birefringence curve, at A = 632. 8nm, T = 
24°C, for: (1) a pure nematic liquid crystal (NLC) mixture, 
denoted as 1294 - 16, of thickness d = 8.08fj.m- and (2) 0.2% 
of antraquinone derivative(AQ) dissolved in 1294 — 16, with 
thickness d = 7.82/imj9j. 
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FIG. 3: Liquid crystal interferometer: a polychromatic plane 
wave, in the ordinary mode, interacts first with a liquid crystal 
grating, with M regions and respective slits, and at the end 
with a lens. A spectrometer measures the intensity of the 
light as a function of the wavelength in the focal plane of the 
lens. 

paths ([3| in the interferometer or the lengths of the M 
paths. The first approach can be implemented by using 
a liquid crystal grating. In the second approach, instead, 
we introduce a generalized symmetric Michelson interfer- 
ometer. 

Let us analyze the first approach. First, we will de- 
scribe an interesting property which makes liquid crys- 
tals able to satisfy the conditions Q for the reproduction 
of truncated exponential sums. When we apply a vari- 
able voltage V to a, liquid crystal cell, interacting with 
an incoming plane wave, we can observe the voltage de- 
pendence of the birefringence An ^ — rio [10' of the 
hquid crystal (see, for example. Fig. [2|. In particular, 
there is a limited range of potentials in which the trans- 
mission, and so the liquid crystal index of refraction, has 
a well defined dependence on the applied voltage. Such 
a definite behavior turns out to be a good tool in order 
to reproduce the terms , with m = 0, 1, M, in Eq. 

The M terms in the truncated exponential sum cor- 
respond, respectively, to M different regions in a liquid 



crystal cell with the same thickness dm = d- 

The basic experimental setup is showed in Fig. |3] 
an incoming polychromatic plane wave, in the ordinary 
mode, interacts with a liquid crystal grating with M 
slits. Such a grating consists of M liquid crystal re- 
gions, with M different variable applied voltages Vm, 
with 171= 1,2, ...,M, and a slit at the end of each re- 
gion. 

So, when the incoming polychromatic plane wave in- 
teracts with the liquid crystal grating, it gives rise to 
M different electromagnetic phase terms, which can be 
manipulated in an appropriate way, varying the applied 
voltages Vm and so the associated indexes of refraction. 
Such terms superpose coherently in the focal point of a 
lens. 

At this point, it is very easy to implement the algo- 
rithm previously described in order to reproduce expo- 
nential sums. 

Of course, in the experimental realization of such ap- 
proach, we need to take into account the dispersion asso- 
ciated with the broadband spectrum of the source. Such 
problem can be overcome performing several measure- 
ments in different ranges of the spectrum of the light 
source such that the relative dispersion in each range is 
negligible. 

It is also important to point out that the larger the 
maximum achievable optical path, the larger is also the 
maximum achievable truncation parameter M in Eq. (|8| . 
Unfortunately, in the liquid crystal approach, the max- 
imum range of variation of the optical paths is limited 
by the thickness d of the liquid crystal cells and by the 
birefringence, calculated when no voltage is applied |llj . 
Consequently, both these parameters determine the max- 
imum number of terms in the truncated exponential sum 
and the maximum range of possible number N we can 
factorize. 

We have seen one possible way of varying the optical 
paths, in an Af-path interferometer, in order to obtain 
phase terms of order j. A second way to achieve the 
same result is varying, in free space, the path lengths dm 
m Eq. In this case, we do not encounter any problem 
associated with dispersion. So we can determine all the 
factors at the same time, exploiting all the spectrum of 
the incoming light [12] . 

This approach can be implemented exploiting the 
multi-path interference in a generalized symmetric 
Michelson interferometer in free space: the usual two- 
path Michelson interferometer is generalized, in a sym- 
metric way, to a M-path interferometer, using M — I 
beam splitters. In Fig. Q we have represented, for sim- 
plicity, the case M = 4 (obviously such an approach can 
be extended to a generic M). In this case, the four inter- 
fering optical paths can be varied arbitrarily by translat- 
ing the mirrors Mi, M2, M3 and M4, respectively. Our 
factorization algorithm can be easily implemented using 
such an interferometer, giving all the factors of any num- 
ber A'^ in a single run. Moreover, because the lengths of 
the interfering paths can be in principle as large as we 
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FIG. 4: Generalized symmetric A/-path Michelson interfer- 
ometer for the realization of exponential sum with truncation 
parameter M = 4. The usual two-paths Michelson inter- 
ferometer is generalized to an M-path interferometer, using 
M — 1 beam splitters. The M = 4 interfering optical paths, 
indicated with dashed, dashed-dotted, continuous, and dotted 
lines can be varied by moving longitudinally the mirrors Mi, 
M2, M3 and M4, respectively, in order to satisfy the condition 



want, there is no limit to the maximum achievable trun- 
cation parameter M and order j of the exponential sum 
we want to reproduce. This is another aspect in favor 
of this approach, rather than the one based on a liquid 
crystal grating. On the other hand such an interferom- 
eter, despite the liquid crystal grating, is more difficult 
to align and to make stable, especially for large values of 
M. Anyway, these challenges can be experimentally over- 
come, and the Michelson interferometer approach turns 
out to be experimentally reliable, especially if we con- 



sider that for exponential sums of high order j, only a 
relatively small number M of interfering paths is neces- 
sary (M - VN). 

In conclusion, we have shown how an optical interfer- 
ometer would allow the factorization of a large number 
A'^ in a single run of a simple algorithm, for any value of 
N. Our algorithm is reliable, in general, for reproducing 
truncated exponential sums of any order j, with a con- 
sequent reduction of experimental resources with respect 
to Gauss sums. We have shown how our approach al- 
lows us to determine, for each value of A^, the discrete 
spectrum of the associated truncated exponential sum, 
as a function of all the possible trial factors. In gen- 
eral, such an approach goes beyond a discrete analysis of 
exponential sums. In fact, new interesting tools in the 
factorization procedure can be provided exploiting the 
physics and mathematics behind the continuous inter- 
ference spectrum of our optical interferometer approach. 
Moreover, we have described two different kinds of opti- 
cal interferometers which can be used in this approach: 
a liquid crystal grating and a generalized Michelson in- 
terferometer. Two separate experiments based on both 
schemes are in progress and the results will be soon avail- 
able in a forthcoming paper. We are also considering the 
use of the same approach for the realization of truncated 
exponential sums with a number of terms M' < A/, ran- 
domly chosen among the total M terms in Eq. [T] 
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